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Abstract: In supersymmetric (SUSY) models with Dirac neutrino masses, a weak-scale 
trilinear Ay term that is not proportional to the small neutrino Yukawa couplings can in- 
duce a sizable mixing between left and right-handed sneutrinos. The lighter sneutrino mass 
eigenstate can hence become the lightest SUSY particle (LSP) and a viable dark matter 
candidate. In particular, it can be an excellent candidate for light dark matter with mass 
below ~10 GeV. Such a light mixed sneutrino LSP has a dramatic effect on SUSY signa- 
tures at the LHC, as charginos decay dominantly into the light sneutrino, vi, plus a charged 
lepton, and neutralinos decay invisibly into z^za We perform a detailed study of the LHC 
potential to resolve the light sneutrino dark matter scenario by means of three represen- 
tative benchmark points with different gluino and squark mass hierarchies. We study in 
particular the determination of the v\ mass from cascade decays involving charginos, using 
the m,T2 variable. Moreover, we address measurements of additional invisible sparticles, in 
our case Z\2, and the question of discrimination against the MSSM. 



Keywords: Supersymmetry Phenomenology, Supersymmetric Standard Model, Dark 
Matter. 



1. Introduction 

There is abundant and compelling evidence that the bulk of matter in the universe is made 
of massive, electrically neutral particles: dark matter (DM) jlj . While the density of DM 
has been precisely determined, the identity of the DM particle (or particles) is a complete 
mystery. Remarkably, many extensions of the Standard Model (SM) of electroweak and 
strong interactions — designed to address theoretical issues related to the breaking of the 
electroweak symmetry — require the introduction of new particles, some of which are excel- 
lent DM candidates. This is most notably the case in R-parity conserving supersymmetry 
(SUSY) g. 

If the origin of DM is indeed new physics beyond the SM, in particular SUSY, there are 
high hopes that DM will be produced in abundance at the LHC, through cascade decays 
of the new matter particles produced in the pp collisions ||. The ambitious goal will 
then be to determine the properties of the DM candidate and reconstruct its thermal relic 
abundance from collider data. On the one hand this may then be used further to make 
testable predictions for direct and indirect DM searches. On the other hand, if the DM 
mass determined at colliders and non-accelerator experiments agrees, it may be used to 
constrain the cosmological model. This would enormously enhance the interplay between 
particle physics, astrophysics and cosmology. The LHC phenomenology of DM- motivated 
SUSY scenarios has hence been discussed in great detail in the literature 0, [| . 

Another experimental evidence that the SM misses something fundamental is the ob- 
servation of neutrino oscillations Q. This can be resolved by including right-handed neu- 
trinos in the model. Current observations, however, do not allow to establish the Majorana 
or Dirac nature of neutrinos. While the smallness of the neutrino mass can be naturally 
explained by introducing Majorana mass terms and making use of the see-saw mechanism, 
Dirac masses for neutrinos with very small Yukawa couplings are a viable and interest- 
ing alternative. In super symmetric theories, one may naturally obtain very light Dirac 
neutrino masses from F-term SUSY breaking ||, In addition to providing an explana- 
tion for neutrino masses, this class of SUSY models offers an interesting alternative to the 
conventional neutralino DM candidate: the sneutrino. 

The crucial point is that in these models one can have a weak-scale trilinear Ay term 
that is not proportional to the small neutrino Yukawa couplings and can hence induce a 
large mixing between left-handed (LH) and right-handed (RH) sneutrinos even though the 
Yukawa couplings are extremely small. The lightest sneutrino can thus become the lightest 
SUSY particle (LSP) and a viable thermal DM candidate. 

Note that the mainly RH sneutrino LSP is not sterile but couples to SM gauge and 
Higgs bosons through the mixing with its LH partner. Sufficient mixing provides efficient 
annihilation so that the mixed sneutrino can be a viable thermal DM candidate with a relic 
density of Qh 2 ~ 0.11 as extracted from cosmological observations ]7|, ||. On the other hand 
the amount of mixing is constrained by limits on the spin-independent scattering cross- 
section, for which the LH sneutrino component receives an important contribution from Z 
exchange; this cross-section is suppressed by the sneutrino mixing angle. Because of the 
gauge and Higgs interactions, the presence of the mixed sneutrino can also significantly 
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impact Higgs and SUSY signatures at the LHC || . 

In H , some of us investigated the case of mixed sneutrinos as thermal DM with special 
emphasis on the mass range below ~10 GeV. We examined the viable parameter space and 
discussed implications for direct and indirect dark matter searches, as well as consequences 
for collider phenomenology. Regarding the latter, we found that the SUSY signatures 
greatly differ from the expectations in the conventional Minimal Supersymmetric Standard 
Model (MSSM) with a Z\ LSP: while squarks and gluinos have the usual cascade decays 
through charginos and neutralinos, with the same branching ratios as in the corresponding 
MSSM case, the charginos and neutralinos decay further into the v\ LSP. In particular, 

%1,2-H'Vi, Wi-tfiui (1.1) 

with practically 100% branching ratio over most of the parameter space. At the LHC, the 
typical cascade decays therefore are q R — > qZ% — > qvv\, q L — > qZ2 — > qvv\ and q L —¥ qW\ — > 
q'l^vi, all giving different amount of missing transverse energy, ftp. Moreover, gluino-pair 
production followed by decays into qq'W\ through either on- or off-shell squarks leads to 
same-sign (SS) and opposite-sign (OS) dilepton events with equal probability. Besides, the 
light Higgs decays invisibly into a pair of LSPs, h° — > v\V\. 

In this paper, we now perform a detailed study of the LHC potential to resolve the 
light sneutrino DM scenario. This includes in particular the determination of the DM mass 
from g — > qq'W\ — > q'l^v\ and/or q — > q'W\ — > q'l^vi events. To this end we rely on the 
subsystem mri method (JuJ. Moreover, we address the question of measuring the masses of 
additional invisible sparticles. The case of a ~ 100 GeV mixed sneutrino LSP was studied 
in 0|. 

The paper is organized as follows. First, in Section [|, we briefly recall the main 
features of the mixed sneutrino model. In Section ^ we present three benchmark points 
and their characteristic signatures. This sets the framework of the analysis. We then go on 
in Section ^| to study the discovery potential at the LHC with 7 TeV center-of-mass energy. 
Measurements at the 14 TeV LHC are discussed in detail in Section [| A summary and 
conclusions are given in Section ^. 

2. Mixed sneutrinos 

The framework for our study is the model of || with only Dirac masses for sneutrinos. In 
this case, the usual MSSM soft-breaking terms are extended by 

A£ soft = m%]Ni\ 2 + A^LiNiHu + h.c. , (2.1) 

where and Ay are weak-scale soft terms, which we assume to be flavor-diagonal. Note 
that the lepton-number violating bilinear term, which appears in case of Majorana neutrino 
masses, is absent. Neglecting the tiny Dirac masses, the 2x2 sneutrino mass matrix for 
one generation is given by 

j J m l + cos W 7l AiV sin 13 \ ra« 
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Here m| is the SU(2) slepton soft term, v 2 = v \ + v\ = (246 GeV) 2 with v\p the Higgs 
vacuum expectation values, and tan/3 = v^/vx- The main feature of this model is that 
the m~, m\ and A„ are all of the order of the weak scale, and Ay does not suffer any 

Li IS 

suppression from Yukawa couplings. In the following, we will always assume < so 
that the lighter mass eigenstate, Ux, is mostly a vr. This is in fact well motivated from 
renormalization group evolution, since for the gauge-singlet Ur the running at 1 loop is 
driven exclusively by the A v term: 

dm 2 - 9 

while 

dm- q q i 

— ^ = -^g\M 2 - -^g 2 Y M 2 + -^A% . (2.4) 

dt 16vr 2y2 2 80rr 2yY 1 16vr 2 " v ' 

A large A v term in the sneutrino mass matrix will induce a significant mixing between 

the RH and LH states, 

. \/2^sin/3 
sm26>j> = — 5 5—, (2.5) 

m~ — m- 

V2 V\ 





leading to mass eigenvalues 



2 



m 2 -„ - f nr. -I- \/ /// 



i + 2,4?t; 2 sin 2 /3^ (2.6) 



where m\ = m|= ztm 2 ^ + m 2 z /2, and < by convention. Notice that a large value of 
A v can induce a large splitting between the two mass eigenstates even if m~ and m 2 ~ are of 

Li JS 

the same order, leading to scenarios where m Ul <C m U2 , . In this way, vx can naturally 
be driven much below the neutralino masses. 

Taking m Ul , m U2 and 9 U as input, the soft terms m^, m^ and A u entering the sneutrino 
mass matrix Eq. ( |2.2D are fixed. This also fixes the corresponding LH charged slepton mass, 
m~ = m~ + m\ cos 2/3(sin 2 9w — 5)- ^ or the RH one, m~ = m 2 ~ — m 2 z cos 2/3 sin 2 9\y, we 
assume = for simplicity. We use an appropriately modified version of SuSpect [12 



for the spectrum calculation, which includes in particular radiative corrections induced by 
the Ay term, as given in 

A full scan of the relevant parameter space was done in S], taking into account con- 
straints from the Z invisible decay width, the Higgs and SUSY mass limits, as well as DM 
constraints from the relic abundance and direct and indirect DM searches. It was found 
that light mixed sneutrino DM consistent with all constraints populates the region 

1 GeV < m ui < 8 GeV and 0.1 < sin 6 V < 0.4 . (2.7) 

Moreover, over most of the parameter space, m U2 > 200 GeV and > ^ . For 

LHC physics this means that the final steps of the cascade decays will be dominated by 
Zi,2 — > v + vx and W\ — > I + v\. The remaining relevant parameters are the squark and 
gluino masses, which determine the production cross sections. 
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SN1 SN2 SN3 
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mu L 
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7U0 


mu R 
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776 
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m ' bl 


1514 


766 


689 


m h 


1441 


675 


584 




811 


807 


805 




524 


524 


524 




516 


516 


516 


m fl 


503 


503 


503 


m Wi 


227 


228 


228 




227 


228 


228 


m Zl 


109 


109 


109 




7.6 


7.6 


7.6 



Table 1: Masses in GeV units for the SNDM benchmark points SN1-SN3. The spectrum is 
computed at the EW scale using a modified version of Suspect |l2|]. For comparison with the 
MSSM case, we consider points MSSM1-MSSM3 with the same masses as SN1-SN3 but the i>\ 
removed from the spectrum. 





7 TeV LHC 


14 TeV LHC 




gg gq qq 


99 gq qq 


SN1 


0.03 0.008 0.0 


1.1 0.5 0.05 


SN2 


0.02 0.2 0.2 


1.0 4.3 2.3 


SN3 


0.002 0.08 0.35 


0.2 2.0 3.2 



Table 2: Cross sections in [pb] for gluino and squark production at 7 TeV and 14 TeV. 



3. Benchmark points and characteristic signatures 

In order to study the light mixed sneutrino DM (SNDM) scenario at the LHC, we pick 
a parameter point from Q with a v\ LSP of 7.6 GeV as the DM candidate. The other 
sneutrino parameters are m^ 2 = 524 GeV, sin 6c = 0.225 and tan/3 = 10 (A u = 348 GeV). 
The neutralino/chargino sector is given by M2 = 2M\ = 221 GeV and [i = 800 GeV. 
Together with = = 514 GeV, this fixes the properties of the weakly interacting 
sparticles. Note we assume flavor degeneracy. 

Based on this setup we define in Table |l] three benchmark points (SN1-SN3) with 
the same electroweak sector but different gluino-squark mass hierarchies. The squark and 



gluino production cross sections, computed using Pythia[|13|, are given in Table and the 
relevant decay branching ratios are given in Table |. 

Employing GUT relations for gaugino-mass parameters, point SN1 has a gluino mass 
of m g = 765 GeV, while squark soft terms are set to 1.5 TeV, resulting in nig ~ 2m g. Point 
SN1 is therefore characterized by heavy quarks and a light gluino, which decays through 
3-body final states into charginos and neutralinos, which will then decay exclusively to X^v\ 
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Table 3: Most important decay channels for points SN1-SN3. 

and vDi, respectively. We hence expect dominantly gluino-pair production leading to 4 jets 
plus 1$t- Moreover, about half of the events will have an isolated charged lepton, and about 
25% of the events will have two leptons with uncorrelated flavor (assuming three roughly 
degenerate light sneutrinos, 2/3 of the leptons are electrons or muons and 1/3 are taus). 
Since the gluino is a Majorana particle, same-sign and opposite-sign dileptons should have 
equal rates. However, same-sign dileptons have less SM background. Promising signatures 
to look for are hence 

(a) 4 jets, leptons, large $t', 

(b) 4 jets, same-sign dileptons, moderate $t'-, 

(c) 4 jets, opposite-sign dileptons, moderate $t- 

Point SN2 has the same gluino mass as SN1 but lighter squarks with niq ~ rrig. 
Therefore SN2 has a much larger overall SUSY production, since the squark-pair and 
gluino-squark associated production give the main contributions to the total cross-section, 
as shown in Table ||. The gluinos have the same decay modes as above, while the squarks 
decay dominantly through q + Wi/Zi t 2- Events with only 2 or 3 jets are predominant and 
result from qq or gq production, with the squark decaying into a neutralino or chargino 
(99% of the ijL and 98% of the q~R decays). 

Finally, point SN3 is characterized by light squarks and heavy gluinos (achieved through 
non- universal gaugino masses). We hence expect dominantly squark-pair production fol- 
lowed by decays into quarks plus charginos or neutralinos, see Table [31 These events have 
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SNlfg^qqZ^ 
SN1 (g -> qq W) 
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MSSM3 (q -» q Z) 
MSSM3 (q -> q W ) 
MSSM3 (q -> q 2) 




Figure 1: Comparison of contributions to the $t spectrum at y/s = 7 TeV for SN1 and MSSM1 
(left) and for SN3 and MSSM3 (right); SN2/MSSM2 gives almost the same picture as SN3/MSSM3. 



2 hard jets plus fir, often accompanied by 1-2 leptons from qr — > qW\ — > qW%. Again, 
events without leptons are expected to have larger Jfir on average than events with leptons. 

It is also interesting to compare the phenomenology of the three SNDM benchmark 
points with the corresponding MSSM cases. To this end we consider points MSSM1- 
MSSM3 with the same masses as SN1-SN3 in Table |l] but with the i>\ removed from the 
spectrum. The MSSM points thus have a neutralino LSP with a mass of 109 GeV. The W\ 
decays exclusively into Z\W, while the Z2 decays into Z\Z (12%) or into Z\h (88% ). 1 The 
larger rate of multilepton events in the SNDM scenario is a clear difference to the MSSM 
case. On the other hand one can expect MSSM events to have higher jet multiplicity. 

Another general feature of the SNDM scenario is that, due to the invisible Z<i decay, 
we expect a harder $r spectrum as compared to a similar MSSM point. This is illustrated 
in Fig. H|, where we show the $ T distribution for the SN1/MSSM1 and SN3/MSSM3 points 
at detector level without any cuts (for details on the event simulation, see next Section). 
Notice moreover that the subset of events with W\ — > + LSP has a slightly harder fir 
spectrum in the SNDM case. This is because the light v\ LSP is more boosted than the Z\ 
LSP in the MSSM case; this difference becomes less significant if the LSP's mother already 
has a large boost. On the other hand, the $t spectrum coming from q/g — > Z\ events are 
identical in both models, since the lightest neutralino decays invisibly in the SNDM. 

Before studying the LHC phenomenology in more detail, a comment is in order con- 
cerning the assumption of flavor degeneracy. As discussed in [||, the exact mass splitting 
between the Pi's of different flavors strongly influences the DM-allowed parameter space. 
There is, however, no difference between one or three light sneutrinos if the mass splitting is 
> 1 GeV. This is because a 1 GeV mass splitting is enough to suppress any co-annihilation 
contributions to Q,^h 2 . Such a mass splitting is easily induced by small splittings in the 
soft terms, which are in fact rather generic even if one starts out with universal soft terms 
at a high scale. If we have, for instance, rriy le > m„ lT , then the heavier v\ e ^ decays 
to the v\ T LSP through 3-body modes [|l4j]. The dominant decay is into neutrinos, while 
visible decays (e.g., v\ e — > e^T v\t) have at most a few percent branching ratio. From the 



1 In the SNDM case, these decays are highly suppressed with respect to the decay into the sneutrino 
LSP. Furthermore, for SN1-SN3, h° — } v\V\ with practically 100% branching ratio. 
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perspective of LHC signatures, the difference to the case of three exactly degenerate light 
sneutrinos is negligible. We therefore take m^ e = = mo lT = for simplicity. Note, 
however, that interesting non-trivial flavor structures beyond the scope of this study may 
appear in the presence of lepton- flavor violation [|l5|, [l6| . 



4. Discovery Potential at LHC7 

We now turn to the discovery potential of the SNDM model at the LHC, with 7 TeV CM 
energy and 0(1) fb _1 of integrated luminosity. The Monte Carlo simulation details are 



presented in Sec. 4.1 and the main distinct signatures are presented in Sec. 



4.1 Event Simulation 

For the SM background, we include in our calculations all relevant 2 — > n processes for 
the multi-lepton and multi-jet searches. Since in this Section we restrict our results to the 
first LHC physics run (< 1 fb _1 and yfs = 7 TeV) we generate (at least) the equivalent of 
1 fb _1 of events for each process, except for our QCD samples (see Table ||). 

For the simulation of the background events, we use AlpGen [lTf to compute the hard 



scattering events and Pythia [13] for the subsequent showering and hadronization. For 
the final states containing multiple jets (namely Z(—t 11, vv) + jets, W(— > lu) + jets, 
bb + jets, ti + jets, Z + bb + jets, Z + ti + jets, W + bb + jets, W + ti + jets and QCD), 
we use the MLM matching algorithm [17| to avoid double counting. All the processes 
included in our analysis are shown in Table |] as well as their total cross-sections and 
number of events generated. The SNDM spectrum was generated with a modified version of 



Suspect [12 1, which includes right-handed sneutrinos, while the decay table was computed 



with CalcHEP p| /micrOMEGAs2 . 4p[ g. Finally, using the SLHA [f(J interface, signal 



events were generated using Pythia6.4 ]lc 

For the event generation, we use a toy detector simulation with calorimeter cell size 
A r/ x A(p = 0.05 x 0.05 and rapidity —5 < rj < 5. The HCAL (hadronic calorimetry) energy 
resolution is taken to be 80%/y/E + 3% for \rj\ < 2.6 and FCAL (forward calorimetry) is 
100%/\/i?+5% for \rj\ > 2.6, where the two terms are combined in quadrature. The ECAL 
(electromagnetic calorimetry) energy resolution is assumed to be 3%/v^E + 0.5%. We use 
the Isajet |2l]] cone-type jet-finding algorithm to group the hadronic final states into jets. 
Jets and isolated leptons are defined as follows: 



Jets are hadronic clusters with \rj\ < 3.0, R = y Ar/ 2 + Ac/) 2 < 0.4 and pr(jet) > 
30 GeV. 

Electrons and muons are considered isolated if they have |r/| < 2.5, pr(l) > 10 GeV 
with visible activity within a cone of AR < 0.2 about the lepton direction, T,p^ lls < 5 
GeV. 

We identify hadronic clusters as 6-jets if they contain a B hadron with px{B) > 
15 GeV, rj(B) < 3 and AR(B,jet) < 0.5. We assume a tagging efficiency of 60% 
and light quark and gluon jets can be mis-tagged as a b-jet with a probability 1/150 
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Cross 


number of 


bM process 


section 


events 


QCD: 2, 3 and 4 jets (40 GeV< Pt(j1) < 100 GeV) 


2.6 x 10 y fb 


26M 


QCD: 2, 3 and 4 jets (100 GeV< Pt{j1) < 200 GeV) 


3.9 x 10 id 


44M 


QCD: 2, 3 and 4 jets (200 GeV< prC/l) < 500 GeV) 


1.6 x 10' fb 


16M 


QCD: 2, 3 and 4 jets (500 GeV< pr(jl) < 3000 GeV) 


9.4 x 10 id 


0.3M 


ti: tt + 0, 1 and 2 jets 


1.6 x 10° id 


5M 


66: oo + 0, 1 and 2 jets 


8.8 x 10 id 


91M 


Z + jets: Z/7(— >■ /t, i/z^J + 0, 1, 2 and 3 jets 


8.6 x 10° fb 


13M 


TXT" i " j_ TT7" - ! - / , 7 \ i r\ 1 O 1 O • J. 

W + jets: W /=c (— >■ tf) + 0, 1, 2 and 3 jets 


1.8 x 10 fb 


19M 


Z + it: Z/7(-> i/P) + + 0, 1 and 2 jets 


53 fb 


0.6M 


Z + 66: Z/j(-^ 11, vv) + bb + 0, 1 and 2 jets 


2.6 x 10 3 fb 


0.3M 


T T T 7 7 T T T— 1— / 7 7 \ i 7 7 i /~\ -1 7 ' "\ ■ i 

W + 66: W /± (-> a//) + 66 + 0, 1 and 2 jets 


6.4 x 10 d fb 


9M 


+ tt: W ± (-^ all) + tt + 0, 1 and 2 jets 


1.8 x 10 2 fb 


9M 


T T 7" i 7 T T T — I- / 7 7 \ "i~7 / 1 7 \ 

W + tb: all) + tb(tb) 


6.8 x 10 2 fb 


0.025M 


tttt 


U.D ID 


11V1 


ttbb 


1.0 x 10 2 fb 


0.2M 


6666 


1.1 x 10 4 fb 


0.07M 


WW: W ± {-+ lu) + W ± (^ Iv) 


3.0 x 10 3 fb 


0.005M 


WZ: W±{-+ lv) + Z{-^ all) 


3.4 x 10 3 fb 


0.009M 


ZZ: Z{-+ all) + Z{-+ all) 


4.0 x 10 3 fb 


0.02M 



Table 4: Background processes included in the discovery potential for LHC7, along with their total 
cross sections and number of generated events. All light (and b) partons in the final state are required 
to have px > 40 GeV. For QCD, we generate the hardest final parton jet in distinct bins to get a 
better statistical representation of hard events. For Wtb production, additional multi-jet production 
is only via the parton shower because the AlpGen calculation including all parton emission matrix 
elements is not yet available. For this process, we apply the cut \m(Wb) — m t \ > 5 GeV to avoid 
double counting events from real ti production. 

for pt < 100 GeV, 1/50 for pt > 250 GeV, with a linear interpolation for 100 GeV 
< p T < 250 GeV. 

4.2 Signal Distributions 

In Fig. ||a-c we show the n(V) and n(j) distributions for the benchmark points along 
with the SM background (BG) after the following set of cuts: 

• $ T > 400 GeV, n(j) > 3, p T (ji) > 150 GeV, p T (j) > 50 GeV and S T > 0.2 

where St is the transverse sphericity and Pt{ji) is the px of the hardest jet. 

Despite having a softer $t spectrum than the SM BG, the point SN1 has a much harder 
n(j) distribution, which peaks at n(j) = 4, as expected from the discussion in Sec. ^. As 
shown in Fig. |2|b, the signal exceeds the SM BG in the n(j) = 6,7 bins. Furthermore, with 
this set of cuts, the lepton number distribution for point SN1 is already at the BG level for 
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Figure 2: $t, n (j), n (l) and ^ Et(1)/ 53 Er(j) distributions for the model points in Table [l] along 
with SM background for LHC7. In frames a ~ c the following cuts have been applied: fir > 400 
GeV, n(j) > 3, pr(ji) > 150 GeV, pr(j) > 50 GeV and S T > 0.2. While frame d has weaker cuts: 
$ T > 300 GeV, n(j) > 2, p T (ji) > 100 GeV, p T (j) > 50 and S T > 0.2 



n(l) = 1, easily surpassing the BG in the dilepton channel. However, due to its small cross 
section, the SN1 signal requires several fb _1 of integrated luminosity in order to become 
visible. We estimate that approximately 5(2) fb _1 are required to claim a 5<t(3ct) evidence 
for the SN1 benchmark. 

As mentioned in Sec. ||, point SN2 has some distinct signatures from the benchmark 
point SN1. While for the latter the sparticle production cross-sections is dominated by gg, 
SN2 has the bulk of its signal coming from qq and qg events. As a result the signal has a 
softer jet distribution, due to squark decays to charginos and neutralinos, as shown by the 
decay branching ratios (BRs) in Table ||. The same is valid for the SN3 signal, which is 
dominated by squark pair production. This is clearly seen in the n(j) distribution shown in 
Fig. |2|b, once the overall signal normalization is taken into account. On the other hand, since 
the 2-body squark decay tends to produce boosted charginos and neutralinos, and because 
of the larger production cross-section, points SN2 and SN3 have a harder ]£t spectrum as 
compared to point SN1 (and the background), as shown in Fig. |2|a. Therefore, the LHC7 
should be able to discover the SN2 and SN3 points with hard $t (^ 300 GeV) and pr(j\) 
(> 250 GeV) cuts in the n(j) = 2 or 3 channels. We estimate that a 5<r evidence for both 
SN2 and SN3 points can be achieved with 1 fb _1 of integrated luminosity. Furthermore, 
from Fig. |2|c we see that both points are also visible in the dilepton channel. 

Since the SNDM signal is rich in boosted leptons, we show in Fig. |2|d the ratio between 
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the scalar sum of the lepton and jet Ex's, Et(1/j) = Yl ^t(0/ S Et(j), for a weaker set 
of cuts: 

• E T > 300 GeV, n(j) > 2, p T (ji) > 100 GeV, p T (j) > 50 GeV and S T > 0.2 

As we can see, both the signal and BG distributions peak at low Et(1/ ' j) values. However 
the SM BG falls much more sharply than the SN2 and SN3 signals, which are above 
background for Ex(l/j) > 0.2. Since the signal and SM distributions have very distinct 
shapes, Ex(l/j) can be used to discriminate between signal and background even for SNDM 
models with smaller signal cross-sections. Nonetheless, point SN1 still is well below the 
SM BG to be seen at LHC7. 

So far each of the SNDM signatures described above are common in standard MSSM 
scenarios, including the CMSSM. Let us now discuss how the distributions shown in Fig. || 
can help to distinguish between the MSSM and the SNDM models. For this purpose we 
also show the respective distributions for a MSSM case (MSSM3) with a spectrum identical 
to point SN3, but without RH sneutrino. 2 As expected, MSSM3 has a softer Et spectrum 
than the SNDM point SN3, although the difference after cuts is not very large. Another 
distinction between the MSSM and SN points is their jet multiplicity. In the SN case all 
charginos/neutralinos decay to v\ + l/v, while this channel is absent for the MSSM points 
and neutralinos and charginos decay instead to the Z\ LSP plus h, Z or W, see Sec. |||. 
Therefore we expect SNDM models to have a softer n(j) distribution when compared to a 
similar MSSM case. This is confirmed in Fig. ^3, which shows that the n{j) distribution for 
the MSSM3 model is suppressed in the n(j) = 3 bin and enhanced in the higher bins when 
compared to the SN3 point (once the total signal normalization is taken into account). 

The main distinct feature of the SNDM scenario appears, nevertheless, in the multi- 
lepton distribution. As discussed in Sec. ||, SNDM models are expected to be rich in hard 
leptons coming almost exclusively from W\ decays. As a result, the leptons in dilepton 
events will have uncorrelated flavors and kinematics. Furthermore, higher lepton multi- 
plicities only appear, at much smaller rates, from top decays. In corresponding MSSM 
scenarios, W\s and Z2 l s decays to leptons are sub-dominant: in case of W\, they are given 
by the BRs of the W, while in case of Z2, they are suppressed by the high BR into h°. 
To give concrete numbers, the ratio of dilepton/0- lepton events is ~ 0.12, 0.06 and 0.06 
for points SN1, SN2 and SN3, respectively, see Fig. |c, while for the MSSM1, MSSM2 and 
MSSM3 points it is ~ 0.04, 0.01 and 0.01. Furthermore, dilepton events in MSSM mod- 
els often come from Z2 decays to l + l~ Z\ and hence consist of same-flavor OS dileptons. 
Therefore the ratio of OS and SS dilepton events is another useful discriminant. For the set 
of cuts used in Fig. 03 we have SS/OS = 0.94, 0.94 and 0.62 for points SN1, SN2 and SN3, 
respectively. The SS/OS ratio is considerably smaller for SN3 than for the other points, 
because its signal is dominated by squark pair production, leading to less SS dileptons. 
Although the SN2 point also has a large qq production cross-section, the n(j) > 3 cut 
suppresses the contribution from qq events to the SN2 signal, enhancing the SS/OS ratio 

2 In order to avoid a proliferation of lines in the plots we do not show the distributions for MSSM1 and 
MSSM2, but we have checked that the general features discussed here by means of SN3 versus MSSM3 hold 
also for the other configurations. 
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in this case. For the corresponding MSSM models and the same set of cuts, we have SS/OS 
= 0.75, 0.54 and 0.39 for MSSM1, MSSM2 and MSSM3, respectively. As we can see, the 
SS/OS ratio is significantly enhanced in SNDM models, when compared to the MSSM. 
Furthermore, as seen in Fig.|2|d, the Et(1/ j) distribution for the MSSM is much softer than 
in the SNDM case, which, together with the n(j) and n(l) distributions can also point to 
a light sneutrino LSP. 

We conclude that, while the Et and n(j) distributions are more model dependent 
and more sensitive to NLO and detector systematics, the n(l) and Et(1/j) distributions 
are promising candidates for an early distinction between the MSSM and SNDM cases. 
Although the first LHC run will probably not accumulate enough luminosity to make use 
of the dilepton invariant mass, depending on the signal cross-section, the dilepton channel 
may already indicate which is the relevant model. However, a more decisive answer will 
have to wait for the later physics run with higher energy. In the next section we discuss how 
LHC14 (y/s = 14 TeV) with an integrated luminosity of 100 fb _1 can give information on 
the SNDM spectrum and thefore provide stronger evidence that the SNDM case is realized. 

5. Measuring Masses at LHC 14 

Despite the distinct features of the SNDM model discussed in the previous sections, in- 
formation on the sparticle spectrum can give a more decisive evidence for sneutrino DM. 
While in most MSSM models, a neutralino LSP consistent with DM and collider constraints 
implies wig > 50 GeV (assuming gaugino mass unification), the SNDM model has a much 
lighter LSP. Since several fb _1 are required for any mass measurement, from now on we 
will present all our results for LHC 14 (^i = 14 TeV and C = 100 ftr 1 ). We include both 
SUSY and SM backgrounds, which include all the relevant processes listed in Table [|, but 
at y/s = 14 TeV (for more details see |22[ ). 

5.1 Method 

Recently much improvement has been made on mass measurement methods at hadron 
colliders and several distinct mass measurement techniques have been suggested, such as 
kinematic endpoints and shapes (invariant mass, s m i n , ...) and transverse mass methods 
(rriT2, TTicr ■>'••)• See |2^, ^4|] for comprehensive reviews and references. Here we adopt the 
subsystem rriT2 method, described in detail in Ref. p0| . It consists in applying the niT2 
endpoint technique to particular subsets of the event topology. 

Gluino-pair production 

First, we consider gluino-pair production, with each gluino going through the cascade 
decay: 

g ->■ qq + W\ -> qq + I + v\ (5.1) 
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From this event topology we can form three subsystems: 

p(l) + P (2) ^ vis (l) + D (l) + vis (2) + D (2) 

i. ~g {l) +g {2) - 



2. + ~g {2) 



(i) 



( 2 ) + Ivf 



(5.2) 



-,(2) 



where the upper index labels the decay branch and the final states are grouped into a 
visible (vis) and a daughter (D) component. For instance, in the first case the gluinos are 
the parents (P), while qql and v\ are the visible and daughter components, respectively. 
On the other hand, for the third subsytem, the charginos are the parents and I and v\ are 
the visible and daughter components. The (subsystem) rriT2 variables are then constructed 
as defined in [25, 26], but using the visible and daughter components defined within each 
subsystem: 



m T 2(m x ) 



mm 



( (1) (2)m 



PT(Dl)+PT(D2) = -PT(visi)-pT (viS2) 



(5.3) 



where 



77? 



(0 



T- y ii>;- i3i +m| + 2( J B T (visi) J E; T (D i )-p T (viSi).pT(D i )) (5.4) 

and m x is the trial daughter mass. Since mT2( m x = wd) < n?p, the value of m^n^mo 
determines the parent's mass in each subsystem. For the above example we have: 
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where and m l T2 are the rriT2 subsystem variables for the subsystems 1, 2 and 



3 defined in Eq. ( |5.2[ ), respectively. Note that for the second subsystem the daughter is 
defined as Wi, which has invariant mass . However, as the above relations show, in 
order to obtain the parent's mass, the daughter mass has to be known. The strength of the 
mr2 method comes from the fact that analytical expressions are known for the function 
m^ x (m x ) and can be used to simultaneously extract both the parent and daughter masses 
from data. For the subsystems 1 and 2 in Eq. (B.2j), we have: 
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The above expressions were derived in Refs. |27| , jlCfl under the assumption of no initial 
state radiation (ISR), so the total parent pt (|pt(<7 ) + Pt(<?^)|) is zero. In Ref. [10] it 
is shown that, unless pt > nig, this is a reasonable approximation. Therefore for m g > 700 
GeV, we can safely neglect the ISR effect. On the other hand, for the subsystem 3, 
the parent system (W^ + W^) is expected to have large px, since they are produced 
from gluino decays and are much lighter than their parents (jn^ = 227 GeV for the 
In this case we have to include the transverse momentum effect in 



cases we consider} 

l,max i 
l T2 



0: 



where ^-t 



_l,max 
T2 



w 1 



(m x ) 



M- + \J (/"- + it) 2 + m l) ~ x > i f ™x < mui 
+ >"i ! - , if rn x > ///,.,. 



Pt \2 



(5.7) 



2m^ 
Wi 
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4m?_ 



Pt 
2m t 



and pt is the chargino pair total transverse 



momentum. Note that, neglecting initial state radiation for the gluino pair, we have 



Pt = \pr(qi) +Pr{q%) + Pr{qz) + Pt(<?4)| = pf 



trans 



(5.8) 

Therefore, selecting events with a fixed p^f ans , Eq. ( |5.7[) can be used to obtain and mo 1 . 
However this cut considerably reduces the statistical significance of the distribution. 

Squark production 

The decay chain (|5.1| ) is relevant for point SN1, where the signal is dominated by gluino 
pair-production, while for points SN2 and SN3, the sparticle production cross-section is 
dominated by squark pair production or squark-gluino production. In this case we will 
consider the following event topology: 



qL ->■ q + Wi -> q + / + v\ 



(5.9) 



The subsystems in this case are analagous to the ones defined in Eq. (|5.2j ), but with one 

less quark in the final state. Using the same notation as before, we label them as m^ 2 , 

q , ; 

ml 



T2 


and m T2 , with: 
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The new mV^ 2 x functions now obey [^7], 10]: 
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= ( + +m? (5-11) 



while m™ ax (m x ) is identical to Eq. (5.7). Note that, for a pair of m x values (m x = and 



m x > mo), the gluino subsystems provide 6 constraints on the 3 masses involved, while 
the squark subsystems provide 5 constraints on 3 masses 3 . Therefore in both cases there is 
an arbitrariness on how to extract the mass values. Since the m l T2 distribution requires an 
extra cut on px in order to allow for the determination of Tn^ and , it will necessarily 
have a much smaller statistical significance than the other m,T2 distributions. Hence we 
will make use only of the first two subsystems, which are already sufficient to determine 
all the masses involved in the process. For that, we adopt the following procedure: 

• For discrete values of m x , we extract the value of m™2 X ( m a;) f° r the subsystems 1 and 
2, from the respective rriT2 distributions, following the algorithm defined in Appendix 

s 

• We then simultaneously fit the results to the appropriate functions, Eqs. (|5.5|), 
Q, O or ( pTD ; 



• The mass values and their uncertainties are then extracted from the best fit result. 

Although only two x measurements at two different m x values are already sufficient 
to obtain all masses, fitting the m™^ (m x ) expressions for a wide range of m x values has 
two main advantages. First, it is less sensitive to uncertainties in extracting m™^ 1 from the 
77JT2 distributions. Second, it allows us to test our underlying model asssumptions since a 
poor fit would indicate that the events selected do not correspond to the topologies (|5j]) 
or (p|). 



Zero-lepton channel 

In principle we can also determine the Z\ and Z 2 masses if we look at the lepton plus 
$T channel. For point SN1 this channel is dominated by gluino decays to neutralinos: 

g-+ qq + Zi/Z 2 ->■ qq + v + h (5.12) 



The m m 2 X (m x ) function for this process obeys Eq. (|5.q ) with — > mj^ 2 . In this case, 
the rriT2 distribution (for a fixed m x ) would present two endpoints, one from g — )■ qq + Z 2 
and one from g — > qq + Z\. Since the latter will necessarily be at higher rriT2 values, a 
simple extraction of w iU always give the Z\ endpoint, with the Z 2 endpoint partially 

obscured by the rriT2 distribution from Z\ events. We have verified that, for the SN1 point, 
the signal/BG ratio is too small in the 4 jets, lepton plus fix channel, making such a 
measurement unviable for C ~ 100 fb _1 . 



3 If, as in Eq. (5.7), we assume a non-zero transverse momentum for the squark pair, Eq. ( |5.1l| ) would 
have two branches and could provide an extra constraint on the masses. However, as in the gluino case, 
ISR effects for squark pair production are negligible. 
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Figure 3: The xtiti distributions for the first two subsystems defined in Eq. (5.2) for point SN1. 
The dashed lines show the signal distributions at parton level, the shaded histogram shows the SM 
contribution and the solid lines show the signal plus background distributions at the detector level, 



after the cuts Eq. (5.14) have been applied. All distributions are normalized to unity. 



For points SN2 and SN3 the lepton plus $t channel is dominated by q decays: 

q — > q + Zi/Z 2 ->■ q + v + v\ (5.13) 



and the m™^) function obeys Eq. ( 5.1l| ) with — > 2 . Since BR(g^ — > Z\) 



100%, while BR(q~L — > Z^) — 33%, the 2 jets, lepton plus channel comes mainly 
from q~R decays and will have a large cross section. Therefore in this case it is possible to 
extract Tn^, once rriq R is known. Hence, after determining the rriq, an( ^ vames 
from q — > W\ + q decays, we will use the zero lepton channel to determine under the 
assumption that left and right-handed squarks are degenerate. 

The results obtained applying the procedure described here to the benchmark points 



SN1-SN3 are presented below in sections |5.2| and |5.3| . In our analysis we include all relevant 
SUSY and SM backgrounds, as well as detector and combinatorics systematics. We always 
assume \/s = 14 TeV and C = 100 fb _1 . 

5.2 Results for SN1 

To apply the niT2 subsystem method for the SN1 point, we consider the event topology 
Eq. ( |5.1| ), which can be selected with the following set of cuts: 

Vt > 100 GeV, n(j) = 4, n(l) = 2, n(b) = 0, 

Mil) > 10° GeV > PrU) > 50 GeV > Pr(0 > 30 GeV . 1 ' ' 

Furthermore, to reduce the Z/j + jets background, we veto OS-dilepton events with 

80 GeV < m(l+r) < 100 GeV or m(/+r ) < 40 GeV . (5.15) 

The total SM BG for this set of cuts is 0.7 fb, while the SNI signal is 4.6 fb. Due to the small 
BG level, our MC samples have only a few events, dominated by ti + 2jets. Although we 
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Figure 4: m™* for point SN1, as a function of the trial daughter mass m x , for the first two 



subsystems defined in Eq. (5.2). The dashed blue lines show the best fit result obtained using 



Eqs. (5.5)-( p.6|) . The solid lines show the exact result obtained at parton level. 



include the full SM background in our results, the effect is subdominant. A more relevant 
issue is the jets and lepton combinatorics. For the subsystem 1 (2), defined in Eq. (p^), 
it is necessary to group the 4 jets and 2 leptons (4 jets) into 2 visible groups, vis^ 1 ' and 
( 2 \ Several methods have been proposed to deal with the combinatorics issue [27, 28]. 



vis 



which usually rely on kinematical correlations between the final states. However, since 
we are only interested in obtaining m™ 1 , for each event we select the grouping which 
gives the minimum value for that event. This way the mn < ni^ x (m x ) relation is 
still preserved even if the wrong grouping is selected. Nevertheless, initial and final state 
radiation (FSR), signal background and detector energy smearing still affects the niT2 
distribution, resulting in a tail for values above m^2 X ( m x)- 

In Fig. ^, we present the an d w|? 2 distributions for the SN1 signal plus the 

SM background, where the trial mass was chosen as the respective daughter mass. For 
comparison purposes we also present the exact parton level distributions (dashed lines) for 
the signal. The spikes in the (solid) rriT2 distributions come from MC fluctuations in the 
SM background. As can be seen from Fig. y, both distributions have an edge at 



• m. 



qql 
T2 



760 GeV and mf 2 ~ 



760 GeV . 



The above values agree well with the expected value, m g . Figure || also shows that the 
and m3p2 distributions are strongly affected by the cuts, ISR, FSR and energy smearing 
effects, when compared to the parton level distributions. Furthermore, our solution to 
the combinatorics usually shifts the m,T2 distribution to lower values, what diminishes the 
peak, resulting in a less evident edge. 

The distributions in Fig. |3| assume a specific value for the trial mass m x , which was 
chosen to be the exact daughter mass in each case. However, as discussed in Sec.||, the 
value of m™2 X as a function of m x allows us to extract both the parent and daughter 



masses, through Eqs. ( p.5|) and (p.6j). Figure |J shows the results obtained from fitting 
the ra™ 2 ai (mi) functions to the values extracted from the simulated data using the 

algorithm outlined in the Appendix. As can be seen from Figs. ||]a and b, the best fit 
for the m|?2 ,maa: and m^" 1 ^ curves (dashed blue lines) agree well with the exact solution 
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777-5 
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QW1 


771 4_ Q 




OQK 4-11 
ZOO ±11 




Q A 4- 1 K 
o4 ± 10 


Exact value 


765 


1520-1523 


227 


109 


7.6 


SN2 




786 ±4 


207 ± 10 


126 ± 13 


14 ±8 


Exact value 


765 


775-779 


228 


109 


7.6 


SN3 




710 ±2 


222 ±5 


111 ±7 


0± 12 


Exact value 


1000 


700-704 


228 


109 


7.6 



Table 5: Measured mass values in GeV for the points SN1-SN3, obtained using the tot2 subsystem 
method described in the text. The error shown only includes statistical uncertainties and assumes 
an integrated luminosity of 100 fb -1 . 



(solid lines). The final result for the three masses, taken as the best simultaneous fit to 
both m™2 X curves, is shown in Table |5|. The approximate precision for rrig and is 
under a few percent, while the precision for the v\ mass is much worse, around 50%, with 
the central value 1.8cr from the true one. This is due to its small mass, when compared 
to the other mass scales, which renders the rri!p£ x expressions weakly dependent on m^. 
Nonetheless, the results clearly point to a very light LSP, with a mass scale much smaller 
than the other particles involved in the cascade decay. 

From the results in Fig. ||| and Table [| we conclude that, despite the lack of precision in 
determining the LSP mass, the 777^2 subsystem method can still show that the LSP state is 
much lighter than expected in most MSSM scenarios. This would provide strong evidence 
for a sneutrino DM scenario, at least for the case of a light gluino/heavy squark spectrum, 
such as in the SN1 point. 

5.3 Results for SN2 and SN3 

For points SN2 and SN3 we must consider the squark cascade decay shown in Eq. (|5.S|). 
To select qi — > W\ + q events we use the following set of cuts: 

y? T > 100 GeV, n(j) = 2, n(l) = 2, 77(6) = 0, 

Mil) > 100 GeV, PT {j) > 30 GeV, p T {l) > 30 GeV 1 ' ' 

and once again veto OS-dilepton events with 80 GeV < m(l + l~) < 100 GeV or m(l + l~) < 
40 GeV. The cross-section after cuts for points SN2 and SN3 are 32 fb and 36 fb, respec- 
tively, while for the SM background we obtain 29 fb. Despite the large background, the 
bulk of the SM events are concentraded at low 777T2 values (< 500 GeV). Thus the position 
of the 777T2 end point extraction is almost unaffected by the SM BG, as seen in Fig. || 
For point SN2 only ~ 65% of the signal comes from squark pair production, due to con- 
tamination from gluino pair production and gluino-squark production. On the other hand, 
the SUSY background is much smaller for point SN3, with ~ 82% of the signal coming 
from squark pair production. After repeating the same procedure used for extracting the 
masses in the SN1 case, but using the appropriate m™^ expressions for the squark decay 
chain, Eqs. ( 5.10| ) and ( 5.11 ), we obtain the best fit result for m^, and shown 



in Table 0. As we can see, for both points the statistical error bars are smaller than the 



-17- 



SN2: m. = 775 GeV, m. = 228 GeV, m. = 7.6 GeV 





m1- 2 (m x = m w ) 




m^Cm = m-) 




SMBG 




Parton Level 



& 0.06 



3 0.04 




550 600 650 
m T2 (GeV) 

SN3: m. = 700 GeV, m. = 228 GeV, m = 7.6 GeV 



n4 2 (m j = m w ) 

mf 2 (m = m^) 

| | SMBG 
1 Parton Level 




550 600 650 700 750 800 850 900 

m T2 ( GeV) 



Figure 5: The m q T2 and va q T2 distributions for the points SN2 and SN3. The dashed lines show the 
signal distributions at parton level, the shaded histogram shows the SM contribution and the solid 



lines show the signal plus background distributions at the detector level, after the cuts Eq. (5.16) 
have been applied. All distributions are normalized to unity. 



ones for point SN1, due to the larger signal cross-section. Once again the least precise 
measurement corresponds to mo 1 , due to its small value. Nonetheless we can still conclude 
that the LSP is much lighter than the chargino and squark. 

As mentioned before, the large qq production cross-section for points SN2 and SN3 
can still provide one more piece of information about the spectrum. From Table t| we have 
BR(q~R — > q + Z\) ~ 100%. Therefore, if instead of q~L pair production we consider the q~Rq~R 
events, we can use the usual rriT2 variable with q~R as the parent, Z\{yv\) as the daughter 
and the jet as the visible component to measure the Z\ mass once rrig R is known. In the 
scenarios considered we have rriq L — rriq R , so we can use the nig value obtained from our 
previous results. To select the right-handed squark signal we require: 

$ T > 200 GeV, n(j) = 2, n(l) = 0, n(b) = 0, p T (ji) > 100 GeV (5.17) 

The SN2 (SN3) signal in this channel is 133 (176) fb, while for the SM BG we have 68 fb. 
Once again, the distribution for the SM background peaks at low rriT2 and has almost no 
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Figure 6: The my 2 distribution for point SN3. The dashed line shows the signal distribution at 
parton level, while the solid lines show the signal plus background distribution at the detector level, 
after the $ T > 200 GeV, n(j) = 2, n(l) = 0, n(b) = 0, p T (ji) > 100 GeV cuts have been applied. 
All distributions are normalized to 1. 



impact on the m™2 x value. In Fig. |6| we show the rriT2 distribution for SN3 with m x = trig , 
where we see a clear edge at rriT2 ~ 690 GeV, very close to the nriq R input value. 
The 

m q,max va i ueg ag a f unc tion of m x for the SN3 point are shown in Fig. [?]. We 
see that the extracted endpoints are very close to their exact value and due to the large 
signal statistics and low BG, the error bars are barely visible. Fitting Eq. ( [5.11 ) to the 
data points in Fig. |?j we obtain: 

2 2 
m QR ~ m 7 

— = 346GeV (5.18) 

which agrees extremely well with the theoretical value, 342 GeV. Assuming rriq L = rriq R 
and using the rriq L value obtained from the qLQL signal (see Table |||), we can compute m% : 

m fi = 111 ±7 GeV (SN3) (5.19) 

where the error in mq R has been included when computing the uncertainty on . Re- 
peating the same procedure for point SN2 we obtain: 



m 



^ = 126 ± 13 GeV (SN2) (5.20) 



The result in this case is worse than for point SN3 because of the larger SUSY background 
present in the SN2 signal. This mainly affects the determination of rriq, as seen in Table [5], 
which propagates to the central value and uncertainty in . Nonetheless, the qnqa 
channel still shows that the measured spectrum has a neutral NLSP with mass ~ m^/2. 
Such a neutral NLSP is consistent with a bino state, if gaugino mass unification is assumed. 
Therefore, the combined results of the qlQl and q~Rq~R channels would point to the usual 
MSSM scenario with gaugino mass unification, but with an additional LSP state, which is 
neutral and very light. This would be another strong evidence for the SNDM model. 
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Figure 7: 



L T2 



for point SN3, as a function of the trial daughter mass m x , in the zero lepton, 



dijct channel, as discussed in the text. The dashed blue line shows the best fit result obtained using 



Eq. (5.11). The solid line shows the exact result obtained at parton level. 



5.4 Dilepton Invariant Mass at LHC14 

From the results in Table |5| we see that the m,T2 method can indicate the presence of a very 
light LSP neutral particle in the signal. Moreover, in case of SN2-3, the m^ 2 distribution 
can give information on an additional invisible sparticle, consistent with the Z\ in case of 
universal gaugino masses. Additional evidence for a sneutrino LSP can be obtained from 



the properties of dilepton events, as already discussed in Sec. [4.2| . Although at 7 TeV 
the dilepton signal is likely too small to allow for the use of the dilepton invariant-mass 
distributions, at LHC14 these may be exploited to probe the nature of the LSP. 

In Fig. H we show the OS and SS dilepton invariant masses for the three SNDM 
benchmark points as well as for the corresponding MSSM models. The cuts applied are 
the same used for the rriT2 analyses, namely Eq. ( 5.14| ) for point SN1 and Eq. ( 5.1(j ) 
for points SN2-3. We assume that the shape of the SM background distributions can be 
extracted from data and/or MC, so we neglect their contribution. As we can see, the mi+i- 
and mi±i± distributions are drastically different between the SNDM and MSSM scenarios. 
Besides the larger overall rate of dilepton events, the SNDM points show a much harder 
mu distribution. This is due to the large W\ — v\ mass gap, resulting in a much harder 
Pt(J) spectrum. Moreover, for the gaugino masses considered here, the Z<i has a significant 
BR to Z + Zx, giving rise to the Z-peak seen in the MSSM OS distributions; in the SNDM 
scenario, no such peak is present since BR(Z2 — > v\ + v) ~ 100%. 



6. Conclusions 

A mainly right-handed sneutrino with large L/R mixing is an excellent candidate for light 
cold dark matter with mass below ~ 10 GeV. If DM is indeed realized in the form of light 
sneutrinos, there are important consequences to the SUSY signatures at the LHC. In par- 
ticular, neutralinos Zx and Z2 appearing in squark and gluino cascades decay invisibly into 
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v\v, so that there can be up to three different invisible sparticles in an event. Charginos, 
on the other hand, decay dominantly into charged leptons plus the v\ LSP. SUSY events 
will therefore present a harder $t distribution than expected in MSSM scenarios with a 
similar sparticle spectrum, and dilepton events will appear at much larger rates both in 
the OS and SS channels. During the first LHC run, at \fs =7 TeV and with ~ 1 fb , a 
signal could already be seen if gluinos and squarks have masses up to ~ 1 TeV. Signal dis- 
tributions such as the lepton and jet number, as well as SS/OS dilepton rates may already 
indicate a light sneutrino as the lightest SUSY particle in the early phase of LHC running. 

Precision measurements enabling a model discrimination should be possible at higher 
energy and luminosity. For y/s =14 TeV and C = 100 fb -1 we have shown that the 
sneutrino mass can be measured using the rriT2 technique with a ~ 50% precision, which 
is already sufficient to distinguish between the SNDM and MSSM scenarios with gaugino 
mass unification. Furthermore, the dilepton invariant mass distributions can also point 
to the presence of a light LSP which carries lepton number. The presence of additional 
invisible sparticles in the decay chains may be inferred from the ftp and transverse-mass 
distributions. We have shown that indeed, for rn^ R ~ m~ L the Z\ mass might be measurable 
with ~ 10% precision. 

Regarding alternative scenarios with possibly similar signatures, a 7-8 GeV Z\ LSP in 
the MSSM with non- universal gaugino masses |2^j (see however |3(J) could be distinguished 
from the case studied here by exploiting, e.g., same-flavor opposite-sign (SFOS) dileptons 
from Z2 — > Z\ + Z, which is absent in the SNDM case. Indeed the absence of kinematical 
structure and flavor correlations typical for the SNDM case will point to W\ — > l^vi 
decays. 4 Quite similar signals as in the SNDM case can in principle arise in the next-to- 
MSSM, with < 10 GeV neutralinos as viable DM candidates that can have large (10 -5 — 
10 -4 pb) elastic scattering cross sections H, H ||. In this case one may have dominantly 
invisible Z2 decays through Z2 — > /i2-^i followed by /12 — > Z\Z\. Here possible ways of 
discrimination are, e.g., the W\ — > W Z\ decays and the presence of additional light Higgs 
states hi and a%. 

We conclude that the LHC offers very good prospects to resolve the light mixed sneu- 
trino DM case. Finally, recall that a corroborating signal is expected in direct dark mat- 
ter searches, as there is a lower limit on the spin-independent scattering cross section of 
a si > 10 -5 pb 

Acknowledgment s 

GB and AL thank the LPSC Grenoble for hospitality. AL would like to thank Xerxes Tata 
for useful discussions, and SK gratefully acknowledges discussions with Sanjay Padhi on 
dilepton searches in CMS. Last but not least, we thank A. Pukhov for providing SHLA 
decay-table output in micrOMEGAs. 

4 In case of a strong hierarchy among the [>i of different flavors there can of course be SFOS dileptons 
from cascades involving Wi'a; however in that case there should also appear the corresponding SFSS events. 



- 21 - 



This research was supported in part by the U.S. Department of Energy, by the Fulbright 
Program and CAPES (Brazilian Federal Agency for Post-Graduate Education), and by the 
French ANR project ToolsDMColl, BLAN07-2-194882. 

A. Endpoint Extraction Algorithm 

Here we describe the algorithm implemented to extract the edge from the mp2 distributions. 
As shown by the solid histograms in Figs.|||, ||| and |(| after the inclusion of both SUSY and 
SM backgrounds as well as ISR, FSR and detector effects, the uiT2 distribution presents a 
tail beyond the expected wS™ x edge. In order to obtain from the mn distribution 

we use a simple linear kink model to fit a subset of the mp2 bins. Although the linear 
kink fit has been used in several rriT2 studies, we have found that in general the result is 
dependent on the region to which the fit is applied. In order to have an unbiased selection 
of the region considered in the fit, we apply the following procedure: 

1. First we select a region consisting of all bins after the peak of the rriT2 distribution, 
such that the initial bin in this region (Mj) is the peak location and the final bin 
(Mp) is the highest mp2 value. 

2. Then we perform a fit to the mp2 distribution in the (Mj,Mp) interval, using the 
linear kink fucntion: 

f , _ J am T2 + b , if m T2 < M 

I [m T 2) - | c( ^ m ^ _ M ) + (1 m + b ,ii mT2 > M y ■ ) 

where a, b, c and M are the free parameters to be fitted. The best fit value for M 
provides an estimative for m^ 1 . 

3. We then select a new region around the best fit M value, consisting of the interval 
(Mj,Mp), where Mi is the same value used for the fit in step 2 minus one bin and 
Mp = 2(M - Mi) + M. 

4. We then repeat steps 2 and 3 until (Mj,M) < 5 bins. 

5. For each of the intervals (Mj,Mp) we compute the chi-square/degrees of freedom 
(x 2 /ndf) for the best fit using Eq. (A.l) and select the (Mi,Mp) interval which gives 



the lowest \ /ndf value. 

6. Finally we fit the mp2 distribution in the selected (Mi,Mp) interval and take mV^^ x 
as the best fit M value. The error on mJp% x is then computed from x 2 — Xmin = 1j 
marginalized over the a, b and c parameters. 

The main idea behind the above procedure is to test distinct mp2 intervals and select the 
one that is better described by Eq. ( A.l| ). Therefore, it provides a general way of selecting 



the proper region to be fitted by the linear kink function. 

Although more sophisticated methods will eventually be used in real data analysis, the 
procedure we have implemented is sufficient for our purposes. Furthermore, we have verified 
that the procedure described here works for different SUSY models and mT2 subsystems, 
despite the distinct shapes of the mT2 distributions. 
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Figure 8: OS (red) and SS (blue) dilepton invariant masses for the SN1-3 points (solid) and the 
corresponding MSSM models (dashed). For frame a) the cuts in Eq.(5.14) have been applied, while 
for frames b) and c) the cuts in Eq.(5.1€) were used instead. 
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